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1. Introduction 

The main objective of the present paper is to revisit our recent contribution and 
give two simpler proofs of the so-called Haff's law for granular gases (with non-necessarily 
constant restitution coefficient). The first proof is based upon the use of Boltzmann's 
entropy and asserts that Haff's law holds whenever the initial datum is of finite entropy. 
The second proof uses only the moments of the solutions and shows that Haff's law holds 



in some weakly inelasticity regime (see Theorem 1.9 for a precise definition in the case of 



non-constant restitution coefficient) for initial datum with finite energy. 

1.1. Motivation. We consider in this paper freely cooling granular gases governed by the 
spatially homogeneous Boltzmann equation 



(1.1) 



\dtf{t,v) =Q,{f,f){t,v) t>0,vGM:' 
\fiO,v) =fo{v), veM.^ 

where the initial datum /o is a nonnegative velocity distribution such that 

/ fo{v)dv = l, I fQ{v)vdv = and / fo{v)\vf dv < oo. (1.2) 

The operator Qe{f,f) is the inelastic Boltzmann collision operator, expressing the effect 
of binary collisions of particles. We assume here that the granular particles are per- 
fectly smooth hard-spheres of mass m = 1. The inelasticity of the collision mechanism 
is characterized by a single scalar parameter known as the coefficient of normal restitution 
^ e ^ 1. Indeed, if v and v^, denote the velocities of two particles before they collide, 
their respective velocities v' and v'^ after collisions are such that 

(u' • n) = —{u ■ n) e, (1.3) 

where n G determines the impact direction, i.e. n stands for the unit vector that points 
from the v-particle center to the t;^-particle center at the instant of impact. Here above 

u = V — Vi, and u' = v' — v'., 
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denote respectively the relative velocity before and after collision. In this work, the resti- 
tution coefficient e is assumed to be a function of the impact velocity, i.e. 

e := e{\u ■ n\). 



In virtue of (O) and the conservation of momentum, the post- collision velocities {v',v'_^) 
are given by 

/ 1+e _ , 1+e _ 
V = V — [u ■ njn, = Vi, -\ — [u ■ n)n. (.1.4J 

The main assumptions on the function e(-) are listed here after (see Proposition |l.ip and 



ensure that the Jacobian of the above transformation ( |1 ■4p is given by 

dv'dvL 



Je{\u ■ n\) 



dvdvi. 



;(|n • n|) + |n • n|e'(|ii • n|) =: 'd'eHu ■ n\) 



where e'(-) and ??g(-) denote the derivative of r h> e(r) and r i-)- 'de{r) respectively (this 
prime symbol should not be confused with the one we have chosen for the post-collisional 
velocity). We refer the reader to [jl]] for more details. The main examples of restitution 
coefficient we shall deal within this paper are the following: 

(1) The first fundamental example is the one of a constant restitution coefficient for 
which e(r) = e G (0, 1] for any r ^ 0. 

(2) The most physically relevant variable restitution coefficient is the one corresponding 
to the so-called viscoelastic hard-spheres [Q] . For such a model, the properties of the 
restitution coefficient have been derived in ^3p and it can be shown that e{z) is 
defined implicitly by the following 

e(r) + ari/^e(r)3/^ = 1 Vr ^ (1.5) 

where a > is a suitable positive constant depending on the material viscosity. 

With the above notations, the Boltzmann collision operator is given, in weak form, by the 
following equation 

Qeifjmv)dv 

\u ■ n\f{v)f{v^) {tp{v') + ip{v'_^) - ip{v) - ip{v^,)) dv dv^, dn (1.6) 

<IR3x§2 

for any smooth test-function i^{v). The strong form of Qe can be recovered easily (see [jl]]). 



1 

2^ 



Notice that an alternative parametrization of the post-collision velocities ( |1.4{ ) would lead 
to a slighty different weak formulation of the collision operator (this alternative formula- 
tion was preferred in [^). 

As explained in [Q] , in absence of any heating source, the granular temperature 

£it)= [ f{t,v)\v\'^dv, t^O 

is continuously decreasing and tending to zero as time goes to infinity, expressing the 
cooling of the granular gases. The precise cooling rate of the temperature is the main 
concern of this note. It was proven in [Q], and predicted by the physics literature long ago, 
that the cooling rate is strongly depending on the choice of the restitution coefficient. Note 
that using the weak form ( |1.6[ ) with the test function %l^{v) = the evolution of £{t) is 
governed by the following relation 



_d 



£{t) = - [ f{t,v)f{t,V,)^e{\u\^)dvdv., (1.7) 
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where the dissipation energy potential associated to e(-) is given by 

r-3/2 /■! 

*g(r) = — - / {l-e^{^/^z)) z^dz r > 0. (1.8) 

We refer to for technical details. In the op. cit. we introduced the following general 
assumptions: 

Assumptions 1.1. The restitution coejficient e(-) is such that the following hold: 

(1) The mapping r G M+ h> e(r) G (0, 1] is absolutely continuous. 

(2) The mapping r G M+ h> i!){r) := r e(r) 15 strictly increasing. 

(3) limsup,,^ooe(?') = cq <1. 

(4) The function x > i — > ^e{x) defined in ( p..8p is strictly increasing and convex over 
(0,+oo). 

These assumptions are fulfilled by the two examples described above and, more generally, 
they hold whenever the restitution coefficient r ^ e(r) is an absolutely continuous and 
non-increasing mapping (see Appendix]). Notice that the first two assumptions are 
exactly those needed in order to compute the Jacobian of the transformation ( ]1.4t ). The 
last two are needed in order to get the following proposition which is based on Jensen's 
inequality. 



Proposition 1.2. Let /o be a nonnegative velocity distribution satisfying (1.2) and let f{t, v) 
he the associated solution to the Cauchy problem Ql.ip where the variable restitution coeffi- 
cient satisfies Assumptions Then, 

and, as a consequence, liuit^oo £{t) = 0. Moreover, if one assumes that there exist a > and 
7^0 such that 

e{z)c^l-az'^ for z ~ (1.9) 
then there exist C > and tQ > such that 

£{t) i^c{i + ty^ yt^to. 



Remark 1.3. The well-posedness of the Cauchy problem ( |1 . ip has been proved in [|Top. 
Remark 1.4. Notice that the above assumption (\1.9\) is equivalent to assume that 



Ue)=sup- — ^<oo (1.10) 

since, for large value o/r > 0, (1 — e(r))/r'^ is clearly finite for any j > 0. In particular, for 
a constant restitution coejficient e(r) = eo, one has £o{e) = 1 — eo which means that dl.lOQ 
holds. For the model of viscoelastic hard-spheres given by ( |1.5p , the restitution coefficient e(-) 
is such that t\i^(e) = a. Furthermore, if £.y{e) < oo for some 7 > 0, then £a{e) = 00 for 
any a / 7. Indeed, the parameter 7 is exactly the one that prescribes the behavior of e{r)for 
small values of r. 



Proposition illustrates the fact that the decay of the temperature is governed by the 
behavior of the restitution coefficient e(r) for small value of r. Now, in order to match the 
precise cooling rate of the temperature and prove the so-called generalized Haff's law, one 
needs to prove that, under Assumptions |1.1| and ( |1.10p , there exists c > such that 
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£{1)-^ c{i + ty^' yt^o. (1.11) 

This was precisely the main objective in [Q] and, as far as the cooling rate is concerned, 
the main result of the op. cit. can be formulate as 
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Theorem 1.5. For any initial distribution velocity /o ^ satisfying the conditions given by 
(|1.2|) with fo G LP{R^) for some 1 < p < oo, the solution f{t,v) to the associated Boltz- 



mann equation ( jl.ip satisfies the generalized Haff's Zaw for variable restitution coefficient 
e(-) fulfilling Assumptions and ( ]1.1C1| ); 

c{i + ty^ i:£{t) s^c{i + ty^, t^o (1.12) 

where c, C are positive constants. 



Remark 1.6. An additional assumption was required in theorem |L5| on the restitution co- 
efficient e(-), see Assumption 4.10]. We do not insist on this point since we believe such 
assumption is only of technical nature and likely unnecessary. 

For constant restitution coefficient, Haff's law has been proved in [|TTt]. This approach was 
generaHzed in leading to Theorem p..5| . The proof is based on the following steps: 



(1) The study of the moments of solutions to the Boltzmann equation using a generaliza- 
tion of the Povzner's lemma developed in [|[] . 



(2) Precise L'p estimates, in the same spirit of OlTO, of the solution to the Boltzmann 
equation for 1 < p < oo. 

(3) A study of the problem ( |l.lp in self-similar variable (that is, for suitable rescaled solu- 
tions) . 

Because of the method of proof, step (2) hereabove, the above result requires strong inte- 
grability assumption on the initial density fo which has to belong to some space with 
p > 1. The main purpose of this paper is to remove the unphysical assumption 

fo € Lf for some p> I, (1.13) 

and prove that the generalized Haff's law still holds under less restrictive assumptions. 

1.2. Main results. We present two independent treatments of the above problem: 



(1) We prove that Haff's law still holds if we replace ( p..l3p by the less restrictive con- 
straint 



and e(-) satisfying ( |1.10| ) for some 7 > 0. 
(2) Using only finiteness of mass and energy on the initial datum, we prove that Haff's 
law ( |4.7p holds true in some weakly inelastic regime defined in the sequel. We notice 



that for constant restitution coefficient, the proof of Haff's law given in [|ri[] required 
the assumption fo G L'p with p> 1. However, it was observed in that, for this case, 
Haff's law holds assuming only that the restitution coefficient is sufficiently close to 
one. We give a complete proof of this fact in the sequel. 

More precisely, the main results of the present paper can be stated in the following theo- 
rems. 

Theorem 1.7. Let e(-) be a non-constant restitution coefficient that satisfies Assumptions 



5| below. Furthermore, assume that the initial distribution fo satisfies ( |1.2[ ) together with 
Ti{fo) < 00, and let f{t, v) be the unique solution to ( |l.lp . Then, the generalized Haff's law 
( |4~7| ) holds true. 

For constant restitution coefficient, our result is weaker, however, we give a qualitative ver- 
sion of Haff's law in this case indicating an algebraic rate of decrease for the temperature. 



see Theorem 3.10 



Referring to the point (2) we state two different results, distinguishing the constant and 
non-constant cases. 
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Theorem 1.8. Let /o be a nonnegative velocity distribution satisfying (1.2) and let f{t,v) 
be the associated solution to ( jl.ip . Assume that the restitution coefficient e(-) is constant 
e(r) = e e (0, 1) and such that 



3(1 



< 1 - «3/2 



where the expression of K3/2 is given in Prop. |2. j[ . Then the generalized Haff's law ( ]4.7D holds 
true. 

For non-constant restitution coefficient, the situation is different and the condition on the 
restitution coefficient will depend on the initial datum. One can formulate our result as 



follows (see Theorem 4.4 for a more precise statement). 



Theorem 1.9. Let /o be a nonnegative velocity distribution satisfying ( |1.2p and let f{t, v) be 
the associated solution to the Cauchy problem ( |l.lp . For any 7 > 0, there exists some explicit 
f-o '■= ^o(/o,7) > such that, if the restitution coefficient e(-) satisfies Assumptions \1.1\ and 
dl.lOD with i^{e) < £0, then the generalized Haff's law ( |4j| ) hold true. 



The proof of Theorem is much simpler than the proof of []2p under the assumption 
( |1.13[ ) on the initial datum. In particular, it does not requires the introduction of self- 
similar variables. It is based essentially on the fact that entropy of the solution f{t,v) to 
dl.lp grows at most logarithmically, namely, there exists Kq > such that 



^(/(t)) ^ifolog(l + t) Vt^O. 

Then, using some estimates which allow to relate the energy £{t) to the entropy, we can 
deduce from such logarithmic growth that the decreasing of the energy £{t) is at most 
algebraic, that is, there exists some finite A > such that inff^o(l + *)^'?(*) > 0. It is 
known from (see also Proposition 2.4) that for non-constant restitution coefficient, this 
is enough to conclude the Haff's law (4.7). The proof of Theorem |1.7| is given in Section 3 



(see Theorem p.9p while several inequalities relating energy and entropy are given in the 
Appendix. 



Concerning Theorems 1.8 and 1.9, their proofs are surprisingly simple and rely only on a 
careful study of the various moments of the solution to the Cauchy problem ( |l.lp . They 
will be the object of Section 4. 



2. Some known results 

We briefly recall some known estimates on the moments of the solution to the Cauchy 
problem. In this section, we will assume that the restitution coefficient e(-) satisfies As- 
sumptions |0| and that the initial datum /o satisfies ( |1.2[ ). We denote then by f{t,v) the 
associated solution to the Cauchy problem (pTTI). For any t ^ and any p ^ 1 we define 



mp{t) 



fit,v)\v\^Pdv 



(2.1) 



with the convention of notation £{t) = mi{t). Then, one has the following proposition, 
see [01. 



Proposition 2.1. For any real p ^ I, one has 
d 

where, 



-mp{t) 



Qe{f,f){t,v)\v\^Pdv ^ -(l-Kp)mp+i/2(t) + Kp 5p(t), 



(2.2) 



[^] 

k=i ^ ^ 
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[2^] denoting the integer part of 2±i and 

is an explicit constant such that Hp < 1 for any p > 1. 

A simple consequence of the above is the following, [Q, Corollary 3.6]: For anyp ^ 1, there 
exists some constant Kp > such that 

mp(0) < oo ^ mp{t) ^ Kp{l + t)~^ \/t ^ 0. 

Furthermore, since we are dealing with hard spheres, the phenomenon of appearance of 
moments occurs in the same way as in the classical elastic Boltzmann [|^ |l^]. Thus, as 
soon as £{0) < oo, the higher moments satisfy sup^^^^ mp{t) < oo for any to > 0. In 
particular, one can rephrase Corollary 3.6]. 



Proposition 2.2. Let fo be a nonnegative velocity distribution satisfying ( |1.2P and let f{t, v) 
be the associated solution to the Cauchy problem ( Jl.lp where the variable restitution coeffi- 
cient satisfies Assumptions |0] and ( |1.9P . For any to > and any p ^ 0, there exists Kp > 
such that 

mp{t) 1^ Kp{l + t)~^ Vi ^ to- (2.3) 
Observe that in order to prove that the second part of Haff's law ( Jl.lip , it is enough to 

3+7 



control 7713+7 (t) in terms oi £{t) 2 . Indeed, recall that 

2 

-£{t)= I f{t,v)f{t,V,)^e{\u\'')dvdv, yt^O. 



_d 



Since ^^(e) < 00, one has 

1 - e(r) ^ ey{e)r^ Vr > 0. 

Plugging this estimate in the definition of *e(''^) and using the fact that 1 — e^(r) < 
2(1 — e(r)) for any r > 0, we get that 



'1 

Jo 4 + 7' 

Since \u\^+^ ^ 2"^+^ {\v\^+^ + \v^\^~^'^), one gets 

-A.£^t) ^ r f^t, v)f{t, V.) {\v\'+^ + d^. 

dt 4 + 7 Jr3xr3 

= ——I^ms_t,{t). 

4 + 7 2 

Therefore, if there exists some constant K > such that 

m3+:^{t) K£{t)^{t) Vt ^ to > 0, (2.5) 

2 

setting C-y = '^^4^^ i-yie), we obtain from ( ^^b that 

—^£{t)!^C^£{t)^ yt^to. 
A simple integration of this inequality yields implies, 

£{t)^ ^ — Vt^to. 

\ + l^£{to)'-^C,it-to)^'^'' 



(2.4) 
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which implies ( Jl.llj ). An additional simplification in the arguments comes with the follow- 
ing proposition which has already been used implicitly in [Q]. We give a complete proof of 
it for the sake of clarity. 

Proposition 2.3. Assume that the restitution coefficient e(-) satisfy Assumptions and is 
such that i-yie) < oo for some 7^0. Assume that there exists some constant C > such that 

m3(t) ^ (t)i(t) VOto>0. (2.6) 
Then, for any p ^3/2, there exists a constant Kp > such that 

mp{t) ^ Kp£{t)P ^ to- (2.7) 

In particular, the generalized Haff's law ( ]4.7| ) holds. 

Proof Let tg > be fixed. First observe that using classical interpolation, it suffices to 
prove the result for any p ^ 3/2 such that 2p G N. Argue by induction assuming that 
for any integer j such that 2j G N, and 1 ^ j ^ p — 1/2 there exists Kj > such that 

mj{t) ^ Kj£{ty for t ^ to. 

Recall that, according to Proposition |2.1| 

d 



^^rnp{t) ^ -(1 - Kp)mp+i/2(t) + Kp Sp{t), 



where 



r p+i 
2 



^p(*) = X] ( ) i^k+i/2{t) rup-kit) +mk{t) mp_fc+i/2(t)) • 
k=i ^ ^ 

The crucial point is that, for p ^ 2, the above expression Sp{t) involves moments of or- 
der less than p — 1/2 except for p = 3/2 which explains its peculiar role. The induction 
hypothesis implies therefore that there exists a constant Cp > such that 

Sp{t)^Cp£{t)p+^/^ yt^to, 

where Cp can be taken as 



Cp = f ^ j {Kk+1/2 Kp-k + Kk Kp_k+i/2) ■ 
k=i ^ ^ 

Furthermore, according to Jensen's inequality mp_^_i/2 

(t) ^ mp+^/^^(t), therefore we obtain 

d 



-^mp{t) ^ -(1 - Kp)ml-^^/^P{t) + Kp Cp£{tY+^'^ yt ^ to. (2.8) 



Additionally, according to ( \3.3\ ) and since e(r) ^ 1 for any r ^ 0, one has clearly *e(|?ip) < 
^ for any n G M'^. Thus, using ( |1.7| ), 

-^-^W ^ ^ / \uffit,v)f{t,v,)dvdv, 

^ / fit,v)\v\^dv = ms/2{t) yt^to. (2.9) 

Let K >he conveniently chosen later and define Up{t) := mp{t)—K£{ty . Then, combining 
(U) and 

±Up{t) = ^mp{t) -pK£{tY-'^^m 

^ _(1 _ Kp)mi+i/2p(t) + Cp £{tr+^/^ + pKmy2{t)£{tr-' Vt ^ to. 
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Therefore, using ( |2.6p we obtain 

This is enough to prove ( |2.7[ ) for K = Kp large enough. Indeed, pick K so that mp{to) < 
K£{to)^. Then, by time-continuity in the moments, the estimate ( |2.7p follows at least for 
some finite subsequent time. Assume that there exists a time > to such that mp{t^,) = 
KE{ti,y, then the above inequality implies 



dUp 



whenever K is large enough. This proves that ( |2.7p holds for anyp ^3/2. □ 

For non-constant restitution coefficient, an interesting result holds: in order to prove the 
lower bound ( p.. lip , it is enough to prove that the cooling of 8{t) is at most algebraic with 
arbitrary rate. More precisely, we have [Q, Theorem 3.7] the following proposition which 
follows from Proposition |2.3| . 

Proposition 2.4. Assume that the restitution coefficient e(-) satisfy Assumptions \1 . 1\ and ( |1.9D 
with 7 > 0. If there exist Cq > and A > such that 

£{t) Coii + t)-^ vt>o, 

2_ 

then there exists C > such that £{t) ^ C {1 + t) i+t for any t ^ 0, i.e. the generalized 
Haff's law ( |4~7| ) holds true. 

3. Entropy-based proof of Haff's law 



The aim of this section is to prove Theorem |1.7| . We begin computing the entropy 
production associated to the Boltzmann equation for granular gases ( p. . IP . 

3.1. Entropy production functional. For any nonnegative /, one can use the weak form 
( |1.6| ) with the test function 'ip{v) = log f{v) to compute the production of entropy 



5e(/):= / Qe(/,/)log/d^;. 

More precisely, 

S,{f) = — I \u. n\f{v)f{v^) log ( ^jjQ^:^) d?;d^;^dn 



2vryiR6x§2 V \f{v)f{v^)j f[v)f{v^) 



+ 77- / k • ^1 {f{v')fiv^) - f{v)fiv*)) (iv dv^ dn. 

^71" -/r6xS2 



Define, 

2vr7M6xs2 V \fiv)f{v^)J f[v)f{v^) J 

(3.1) 

which is a non-negative quantity since log x ^ x — 1 for any x > 0. Notice that, if e = 1, 
then De is the classical entropy production functional and Se{f) = —Ve{f) ^ 0, which 
means that the entropy production Se is non-positive. 

For inelastic collisions, the entropy production functional is more intricate, 

S,{f) = -Ve{f) + ^ I \u-n\{f{v')f{v'J-f{v)f{v.)) dvdv.dn 
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which means that the entropy production splits into a dissipative part (— Pe) and a non- 
negative part. Let us compute more precisely this last term. Since i?e(-) is strictly increas- 
ing, it is bijective. Moreover, \u' ■n\ = -i^edu • n|)), thus, one can write \u-n\ = "l}~^{\u' ■n\). 
Then, using the change of variables (w', v'^) — ;> {v, v^) we obtain 

\u-n\f{v')f{v'Jdvdv^dn= [ ^-^{\u' ■ n\)f{v')f{v'^) dv dv^dn 

j-i.i N ^ dvdvi,dn 



^-\\u-n\)fiv)f{v.) , . 
<§2 Jg(t?e 

We used that 

dv' dv'^ = Je{\u ■ n\) dv dv^, = Je(i?~^(|n' • n|) dv dw*. 
It is easy to see that t?r^(|n • n\) = , IT"^ . Then, we deduce that 

— /" In • n| f (i;) f (i;^) f 1| dudw^^dn 



VAf{v)f{v^)^e^u\)dv dv^. 



For any fixed f j,, we have defined 



u -nl \ 1 1 1 I dn 



ei^7Wu-n\))Jei^7'{\u-n\)) 



After some minor computations. 



Setting then r = "de^i^) ^i^d recalling that Je{y) = i^'eiu): we easily get that 

*e(kl) = ^ / (r- - i?e(r) <(r)) dr 

where we also used that i^e^i^) = 0. We just proved the following proposition. 

Proposition 3.1. Assume that the restitution coefficient e(-) satisfies Assumption \1 . 1\ items 
(1) and (2). Then, for any non-negative distribution function f{v) 

Se{f)= f Qe{fJ)\ogfdv = -Ve{f)+f \u\f{v)f{v,)^e{\u\)dvdv, (3.2) 



where Ve{f) ^ is given by ( p^ while *e(-) is defined by 

^e{g) = ^ {r-Mr)K(.r)) dr, ^gyO. (3.3) 

Q'' Jo 

Additional qualitative properties of $e are given in the following lemma. 

Lemma 3.2. Assume that the restitution coefficient e(-) satisfies Assumption \1.1\ Assume 
moreover that ( |1.9[ ) is satisfied for some a > and 7 > and that there exist two positive 
constants C > and m ^ 1 such that 

i9^Hy) ^ Cy"" for large y. (3.4) 

Then, ^ed""!) ^ C\u\'^^'^^^^ for large \u\, and $e(|u|) ~ 2a\u\"' for small \u\ ~ 0. 
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Proof. Since t?e(-) is assumed to be increasing, one clearly has r — de{r)'d'^{r) ^ r for any 
r ^ 0. Therefore, 



and the first part of the Lemma follows from ( [3.4{ ). Moreover, if e(r) ~ 1 — ar'^ for r ~ 
and 7 > 0, then r — ^eir) "d'eir) ~ a{2 + ^)r"'+^ for r ~ 0. Since 'i?^"'^(r) ~ r for small r, one 
gets easily the second part of the result. □ 

The case of a constant restitution coefficient is included in the previous lemma, however. 



in this case $e is explicit, we refer to [ ]TT| , |^] for previous uses of the entropy production 
functional in the constant case. 

Example 3.3 (Constant restitution coefficient). If e{z) = e G (0, 1] for any 2^0, then 
Je = e and 

^ , , 2(l-e2) f^/- l-e2 

^e{g) = 2 / rdr = 2— • 

Q'^ Jo 

Example 3.4. If e{-) is the restitution coefficient for visco-elastic hard-spheres, there exists 
a> such that 

1 3 

e(r) + ar5 e5 (r) = 1 Vr > 0. 

5 _i 

One checks without difficulty that e(r) ~ a 3r 3 as r oo. /n particular, 

_5 2 

t^eCr) ~ a 3r3 as r oo. 

3 

Consequently, there exists some positive constant Ca > such that ^(y) ^ Coy 2 for large 



y > 0. Therefore, the assumption (3.4) of the above Lemma is fulfilled with m = 3/2, and one 



obtains that there exists some positive constant C > such that 

^ei\u\) < ^1^1 for large u £ M^. 



Since ( |1.9p is known to hold with 7 = 1/5 and a = a, we also have 

*e(|n|) ~ 2a|u|^/^ as ~ 0. 

3.2. Evolution of the entropy and the temperature: Haff's law. In all this section, we 
shall assume the following additional conditions on the restitution coefficient. 

Assumptions 3.5. Assume that the restitution coejficient fulfills Assumptions Moreover, 



assume that (1.9) and (3.4) holds, that is. 



(1 ) There exist a > and 7^0 such that 1 — e(r) ~ ar^ as r ~ 0. 

(2) There exist m ^ 1 + 7/2 and C > such that i^e^iu) ^ Cy"^ for large y. 

Remark 3.6. Note that the assumption m ^ 1 + is no restrictive since the condition (2) 
concerns large values ofy. 

Under this conditions, the growth of the entropy of the solution to ( |l.lp is at most loga- 
rithmic. 

Proposition 3.7. Assume that the restitution coefficient e(-) satisfies Assumptions |3.5[ In 
addition, assume that the initial distribution /o satisfies ( \l.2\ ) together with 'H(/o) < 00 and 
let f{t, v) be the solution to Then, there exists a constant Co > such that the entropy 

n{f{t))of f{t,v) satisfies 

n{f{t))^n{fo) + Coiog{i + t) yt^o. 
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Proof. From the results of previous section, the evolution of the entropy V.{f{t)) is gov- 
erned by 

4-Mf{t)) = -V,{f{t))+ f \u\f{t,v)f{t,v,)^e{\u\)dvdv, Vt^O. (3.5) 

Under Assumption p.5[ , Lemma |3!^ implies that we have *e(|n|) ^ C|np*^™~^) for large \u\ 
while $e(|^i|) — 2a|u|''' for small \u\ ~ 0. In particular, there are two positive constants A 
and B such that 

*e(|u|) ^ ^|n|T + B|n|2(™-i) V-uGM^ 
From dsTsl ) and since -Ve{f{t)) ^ 0, 



Moreover, since ^ 2'^ + 1^^*^^^) one has 

H'^'^'^fit, v)f{t, v^) dv dv^ ^ 2^+1 / br+V(i, = (t) 



Moreover, setting p = m — 1/2, one notices that there is a constant Cp depending only on p 
such that 

/ |up™-i/(t, v)f{t, V,) dv dv. ^ mp{t) 



where the mp{t) terms are the p'^^ order moments defined in ( ]2.ip . Using Proposition |T]2 



together with Proposition |2.2| , one concludes that there exist two positive constants Ci and 
C2 such that 

^^(t) ^Ci(l + t)-i + C2(l + t)"^ Vt>0. 
Since p = m — 1/2 with m ^ 1 + 7/2 one has ^ 1 and, setting Cq = Ci + C2, we get 

-^n{t) Cq{i + 1)-^ vt>o 

which yields the conclusion. □ 

The above logarithmic growth is exactly what we need to prove Haff's law. Indeed, the 
following general result allows to control from below the temperature using the entropy. 
The proof of the following proposition is given in the Appendix. 

Proposition 3.8. Let C denote the class of nonnegative velocity distributions f = f{v) with 
unit mass, finite energy and finite entropy 



f{v)dv = l / Itj] /(u)dv < 00, T-L{f)= / f {v)\og f {v) dv < 00 



Define 



H(/)= / f{v)\\ogf{v)\dv, fee. 
Then, there is some constant c > such that 



n{f)^n{f) + c( [ 



f{v) \v\ dv 



\ 5/3 



and some other constant C > such that 



f f{v)\v\^dv^Cexp{-lH{f)) V/gC. (3.6) 

JR3 



Proposition |3.8| combined with the logarithmic growth of T-L{f{t)) prove the Haff's law for 
non-constant restitution coefficient. 
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Theorem 3.9. Assume that the restitution coejficient e(-) satisfies Assumptions ^.5\ with 7 > 
0. Let the initial distribution fo satisfies ( |1.2p and H{fo) < 00, and let f{t, v) be the unique 
solution to ( |l.ip . Then, there is a constant c > such that 

In particular, the generalized Haff's law ( ]4.7D holds true. 



Proof With the notations of the above Proposition p.8| , there is some constant C > 
independent of time such that 

Since £{t) is bounded, one can find constants ko,Ko > such that 

H(/(t)) ^ ko + Ko log(l + t) Vt ^ 0. (3.7) 



Proposition p.8| also impHes that there exists a constant ci > such that 

£{t) ^ ci exp (^-lH(/(t))^ ^ C2 exp (^-^ log(l + t)) Wt ^ 0, 

with C2 = ci exp(-^). Therefore, setting Aq = ^ > we get that 

£{t) ^ C2(l + 1)"^^" Vt ^ 0. (3.8) 



Therefore, according to Proposition |2.4| , the estimate ( p.SP is enough to prove the second 
part of Haff's law (jTTlp . □ 

For constant restitution coefficient, our result is less precise, however, we prove an inte- 
grated version of Haff's law in the next theorem. 

Theorem 3.10. Let e G (0, 1) be a constant restitution coejficient and let the initial distribu- 
tion fo satisfies ( |1.2| ) with 'H(/o) < 00 , and let f{t, v) be the unique solution to ( |l.ip . Then, 
there are two positive constants a,b > such that 

ft , 1 

/ y^ds ^ -log(l + a6t) Vt^O. (3.9) 
Jo " 

Consequently, 



sup<^ A ^ 0, sup{l + ty^£{t) < 00 ^ = inf <^ A > 0; lim sup{l + ty^ £{t) > ^ = 1. 

Proof Recall that, for constant restitution coefficient e G (0, 1) the evolution of the entropy 
is given by 

4-Mfit)) = -'De{f{t)) + ^—^[ \u\f{t,v)f(t,v,) dvdv. Vt^O 

where we used Eq. ( |3.5[ ) and Example ^]3|. Arguing as in the previous proof, we see that 

-fn{f{t)) ^ 2i^mi/2(t) ^ 2^-2^ Vt ^ 0. 
Integrating this inequality yields 

n{f{t)) ^ n{fo) + r Vrnds yt ^ o. 

e Jo 

Consequently, there exists a positive constant A'l > such that 

H(/(t)) ^ ETi + / y^ds Vt ^ 0. 

e JO 
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Then, from Proposition p.8|, 



£{t) ^ Cexp (^-^ - ^^^3/'^ 1^ V£{s)ds^ Vt^O, 
i.e. there are two positive constants a = \fC exp(— ^y^) and h = such that 



^/£(t) ^ aexp (^-b ./liTjds^ Vt ^ 
which yields (|3.9p . Now, setting 



^ = {A ^ ; sup(l + tf^£{t) < 00} 



we know from Proposition p. .21 that 1 G A, i.e. A ^ 0. Then, it follows from ( |3.9D that 
sup^ = 1. Now, let us define 



B = {X^O; limsup{l + ty^£{t) > 0}. 



Notice that inequality ( |3.8[ ) holds for constant restitution coefficient. In particular, it proves 
that B ^ 0. Notice also that if Ai G B, then any A2 ^ Ai belongs to B. We argue 
by contradiction to prove that inf B = 1. Otherwise, from the previous observation, one 
would have inf 5 = Ag > 1. Pick A € (1, Ag), it follows that X ^ B, that is, limsup^^oo(l + 
tY^8{t) = and, in particular, \ G A. This is impossible since sup^ = 1, and thus, 
infe=l. □ 

Remark 3.11. Though less precise that the converse inequality ( Jl . 1 ip , the above integrated 
version ofHaff's law asserts that (1 + t)^^ is the only possible algebraic rate for the cooling 
of the temperature £{t). Notice that vrit) = \/£{t) is proportional to the so-called thermal 
velocity ^ and we may wonder what is the physical relevance of the above identity ( ]3.9[ ). 
Finally, we recall that it is expected the existence of a self-similar profile ^h{-), homo- 
geneous cooling state, such that f{t,v) = v-r(t)~^$f/( ^ ^.^x ) is a solution to The 
existence of homogeneous cooling state has been proven in (with a slightly different def- 
inition where vr(t) was replaced by (1 + t)^^) and where the self-similar profile is, by 
construction, satisfying £ for some p > 1. We conjecture that the existence of such an 
homogeneous cooling state can be obtained using only entropy estimates, that is 1-L{(^h) < <x>. 

2_ 

Remark 3.12. Notice that, since £{t) ^ C(l + t) 1+7, it follows from ( |33[ ) that there is 
some constant K > such that 

U{f{t)) ^ Klog{l + t) Vi^O 

which means that the logarithmic growth obtained in Proposition |3.7| is optimal. 



4. Haff's law in the weakly inelastic regime 



This section is devoted to the proof of Theorems |T]§ and [L^- Let us now explain briefly 
the strategy of proof to get a precise version of generalized Haff's law. Note that due to 
Proposition |1.2| , one only has to prove a lower bound of the type 

£{t) ^ c{l+t)~^ 

for some positive constant c > independent of time. The following approach uses only 
the evolution of some moments of the solution f{t, v) with the particular use of Proposition 



2.3. We will distinguish between the case of a constant restitution coefficient 7 = and 



the non-constant case 7 > since the two results are different. 
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4.1. The case of a constant restitution coefficient. We assume here that the restitution 
coefficient e(-) is constant: e(r) = e for for any r G M+. In this case, 

1 - / 

*g(r) = Vr > 

8 

and the evolution of the temperature is given by 

4^(t) = -i^/" f{t,v)fit,v,)\v-v^fdvdv^, t^Q. 

Since \v — fj.^ ^ {\v\ + \vi,\f' = \v\^ + + 3|?jp|z;*| + 3|w||t)^p, one deduces that 

- ^/{t) ^ i^m3/2(t) + m^i/2{t) yt ^ 0. (4.1) 

This observation will help us in proving the following theorem. 

Theorem 4.1. Let fo be a nonnegative velocity distribution satisfying ( |1.2p and let f{t, v) be 
the associated solution to ( Jl.lp . Assume that the constant restitution coefficient e is such that 

fc^<l-.3/2. (4.2) 
Then, for any to > 0, there is an explicit positive constant Cq > such that 

m3/2{t) i^CoSitf^^ yt^to. (4.3) 
Consequently, the second part ofHajf's law holds 

£{t) ^ ^ 2 ^ 0. 

\ + Co ^/m{'^ - e2)i' 

Proof Let io > be fixed. According to Proposition |2.1| 

-^"l3/2(t) ^ -(1 - «;3/2)m2(t) + K3/2 S'3/2(0- 

From the expression of 53/2(0 one gets 

-^m3/2(t) ^ -(1 - K3/2)"i2(t) + "l3/2(t)mi/2(t) + £'^{t) Vt ^ to (4.4) 

where we used the fact that K3/2 < 1- Let K be a positive number to be chosen later and 
define 

Uz/2{t) ■.= m^/2{t)-K£{tfl\ 

From Holder's inequality, 

"^3/2 (t) ^ \/W)^m2{t) and mi/2(t) ^ Vt ^ 0, (4.5) 

so that 



Now, using (|4.lD, one gets 

nlJt) 



df7o/2 "^3/2(0 I o 



dt ' ' ' £{t) 



+ " ^ Km^/2{t)£{tf/^ + ^^\''^ K£itf. yt > to. (4.6) 
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This last inequality, together with the smallness assumption ( |4.2p imply the result provided 
K is large enough. Indeed, choose K so that m^/2{to) < K£^/'^{to). Then, by time- 
continuity of the moments, the result follows at least for some finite time. Assume that 
there exists a time > to such that m3/2(t,t) = K£'^/'^{U), then ( |4.6p implies 

where a := (1 — K3/2) — |(1 — e^) > by (|43|). Choosing K = Cq large enough such that 
'^^^/'^ (t^) ^ for any the conclusion holds. □ 



dt 



Remark 4.2. We just proved that, under the sole assumption ( [4.2p , the solution f(t, v) to the 
Boltzmann equation ( p..lp satisfies the generalized Hajf's law 

c(i + t)"2 ^ ^ c(l + ^)"^ t^o (4.7) 

where c, C are positive constants depending only on £{0). 

Remark 4.3. L/sing the explicit expression 0/^3/2 given by Proposition [2. j[ one sees that ( |4.2[ ) 
amounts toe^ \/^^ - f - 0.809. 

4.2. Non-constant case 7 > 0. We consider in this section the non-constant case 7 > 
with restitution coefficient e(-) satisfying Assumption p. . 1| with 

1 - e(r) 

i'y(e) = sup < 00 for some 7 > 0. 

r rT 



The proof is more involved but still based on Proposition |2.3| . We therefore need only to 
estimate m^/2{t)- We will work with the following class of initial datum: Let Sq, qq be two 
fixed positive constants. Define ,^{£o,qo) as the set of nonnegative distributions g ^ L\ 
such that 

/ g{v)<lv = l^ I g{v)vdv = {), I \v\^ g{v) dv = £q 

Jrs Jr3 Jk3 

and 

3/2 

g{v)\v\^ dv Qo ( I g{v)\v\'^dv 



With this definitions we have the following result. 



Theorem 4.4. Let /q be a nonnegative velocity distribution satisfying ( |1.2p with initial energy 
£0. Then, there exists £0 := £o{£o,j) > such that, if the restitution coefficient satisfies 
£-y{e) < £0, then solution f{t, v) to ( |TTT| ) fulfills the generalized Haff's law <^%. 

Proof Throughout the proof, we simply denote ^^(e) = A since both 7 and e(-) are fixed. 
Fix a time to > 0, due to appearance of moments, there exists a e-independent constant 

C := C{tQ, £q,^) > such that sup^^t^ ^2^(1) ^ C. Note that /(to, v) € ^{£tQ, Qto) where 

"13/2(^0) 

— — < Qto < 00. 

mi(to)^/^ 

Additionally, recall that 

-^rn3/2(t) ^ -2am2(t) + m3/2(t)mi/2(t) + £^it) Vt ^ to, (4.8) 
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where we have set 2a = (1 — K3/2) > 0. Now, using (2.4) 

dt 



For t ^ to the quantity m2-Y(t)^/'^ is controlled uniformly by A := C^/'^, therefore, 

^ (f(t)2 + 3^m2(t)) Vt^to- (4.9) 

Set C/3/2(t) := "13/2(0 ~ K£{t)^/'^ where the constant K > will be suitable determined 
later on. Use estimates (Ot) and (K% to get for all t ^ 



d 3c 
^f^3/2(i) ^ -2Qm2(t) + m3/2(t)mi/2(t) + -^XK£{tf 

+ ^XKAm2{t) 
o 

^ -am2(t) + m3/2(t)mi/2(t) + (^1 + Ai^^^ £{tf 

+ f-a + Ai^^^^ m2(t). 



For the first term in the left-hand side, one uses the same estimate as in Theorem 

-am2(t) +m3/2(t)mi/2(t) ^ -a£{ty^ml^{t) + my2{t)£{t)^^\ 

thus, 

d 

~dt 



-U^,2{t) ^ -a£{t)-^ml/^{t) + m^/2{t)£{tf'^ + (1 + 3Kc^) £{tf 



+ {-a + XK^A\m2{t) Vt ^ to- (4.10) 



Without loss of generality, we have assumed A < 8 (any other positive number would have 
worked). Now, let i^o > be the positive root of -aX"^ + (1 + 3c^)X + 1 = 0. Then, if 
K = max(Ko, Qt^^), one gets that 

A < 4 ^ m^/2{t) ^ K£^/'^{t) Vt ^ to 

where = 9^^° k ■ Indeed, since C/3/2(to) < 0, by a continuity argument, then U^/2{t) 
remains nonpositive at least for some finite time. Assume that there exists a time t* > to 
such that U^/2{U) = 0, that is, 7713/2 (t^) = K£^/'^{t^). Then, from ( H.lOp 

^f^3/2(i*) ^ {-aK^ + {l + 3c^)K + l)£{t,f+(^-a + XK^A^ m2(t,). 

Since K ^ Kq, the first term on the right-hand side is negative while, by choice of A < 
^0 := io{to,£o,j), the last term is negative as well. In other words, -^C/3/2(i*) < from 
which we deduce that U^/2{t) will remain non-positive for all t ^ to. As explained in the 
paragraph precedent to Proposition |2.3| , this yields Theorem |4.4] . □ 

Remark 4.5. It is important to notice that the smallness condition is depending on the initial 
datum /o- This is the major difference with respect to the constant case where the smallness 
assumption ( |4.2[ ) is universal. 
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Appendix: Functional inequalities relating moments and entropy 

In this section, we present some functional inequalities that relate moments and entropy. 
We present the results in M" with n ^ 1 regardless we only use them in dimension n = 3. 
In this section, for any measurable subset E c M", \E\ will stand for the Lebesgue measure 
of E. Let / = f{v) be a nonnegative distribution and denote 



nif)= fiv)\ogfiv)dv, H(/)=/ /(^;)|log/(^;)|dT;, 
while, for any k ^ we set 

Mfc(/)= / f{v)\v\Uv. 
First recall the following simple estimate which can be traced back to []7p : 

Lemma A. 1. If there is k £ (0,n) such that Mf,{f ) < oo and < oo then H(/) < oo. 
More precisely, for any k G (0, n) there exists some positive constant Cn^k > such that 

n 

O^U{f)^n{f) + Cn,kMk{f)n+k, (A.1) 

Proof Although it is almost explicitly stated in [^, we provide the complete proof of the 
above estimate. Define A = {v gW\ f{v) < 1}, thus 



H(/) = / f{v) log f{v) dv - / f{v) log f{v) dv = n{f) - 2 / f{v) log f{v) dv 
JR"\A Ja J a 

= «(/) + 2 /^/W log (^)d. 

For any a > 0, let B = e M" ; f{v) ^ exp(-a|v|'')} and B'' = W \ v e Ar\ B then 
log(j^) ^ a\v\^ and 

H(/) ^ nU) + 2a [ fivM" dv + 2 [ f{v) log dv 
Jahb Js-nA \f{'")J 

^ nif) + 2aA4(/) + 2 / f{v) log f-i 



dv. 

Since xlog(l/x) ^ for any x G (0, 1), where M = 2exp(— 1), we get 



f(y) log d^; ^ M^^^ ^/W)^v < M J^^ exp (-«^) di-. 

Set 

n 

where r(-) is the Gamma function. We have H(/) ^ Hif) + 2aMk{f) + 2J„,fc(a) for any 
a > 0. Therefore, for any k G (0, n), there exists some positive constant C„ ^ > such that 

H(/) ^ nU) + 2aMk + Cn,ka'^ Va > 0. 
Optimizing the parameter a > yields ( |A..lp for some explicit constant c„ fc. □ 

We now give a general estimate which allow to control A/fc(/) from below in terms of T-L{f). 
It is likely that such an estimate is well-known by specialists. We include a proof below. 

Proposition A. 1. Let / ^ be such that H(/) < oo with fiv)dv = 1. Then, for any 
k ^ and any e > 0, there exists C{n, k,e) > independent of f such that 

Mkif) ^ C{n,k,e)exp (^---^A_H(/)^ . (A.2) 
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Proof. For any R> 0, let Br denote the ball with center in the origin and radius R, and let 
B'^ be its complement. Then, 



Mfc(/)= / f{v)\vfdv+ / f{v)\v\'dv 
Jbr Jb'j^ 



R'' I f{v) dv = R''(l- [ f{v) dv^ . (A.3) 



Therefore, in order to control Mk{f) from below it is enough to control the mass of / on a 
suitable ball. Let us fix i? > and recall the generalized Young's inequality 

xy ^ X log X — X + exp(y) Vx > 0, Vy G M 

or, equivalently 

xy ^ - log - - - + exp(A y) Vx > 0, A > 0, My £ R. 
A A A 

Using this inequality with x = f{v) and y = I and integrating the inequality over Bji we 
obtain 

[ f{v)dv^^ [ f{v)logf{v)dv-UlogX + l) [ f{v)dv + expiX)\BR\ VA > 0. 

J Br ^ J Br ^ J Br 

Consequently since J^^ f{v) log f{v) dv ^ H(/), we get 

i(logA + A + l) /" f{v)dv^\^{f) + eMX)\BR\ VA > 0, 
^ Jbr ^ 

that is, 

(log(Aexp(A)) + l) / f{v)dv^n{f) + \eMX)\BR\ VA > 0. 
Jbr 

Set X = log(Aexp(A)) and assume x > —1. The last inequality becomes 

^ H(/) + exp(x)|i?fl| 
f(v)dv!^ Vx > -1. 

Br '' 1+^ 

We optimize the parameter x > — 1 by noticing that the right-hand side reaches its minimal 
value for xo such that 

xoexp(xo)|Sij| = H(/). 

The mapping x i-)- xexp(x) is strictly increasing over (—1, oo) and we define W its inverse 
(Lambert function). We get then xq = W (H(/)/|5/j|) and 

f ^ H(/) + exp(xo)|i3fl| f,,,f^(f)\\ir,\ 

Combining this with (|A.3p we get 



Mfc^i?^(^l-exp(^Ty(^^jj|B^|j Vi?>0 

and we still have to optimize the parameter R. For simplicity, set Y = H(/) and X = \Br\. 
For any e G (0, 1), 

1 - expiW{Y/X))X = £ ^ WiY/X) = log (^] ^ Y/X = ^ log " ^ 
where we used that W^^{t) = texp(t). Thus, 



(Him 



X ' X °\ X 



1 - exp(Ty(y/X))X = e ^ ^ = log (^) ^ X = (1 - e) exp ^ ^ 
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To summarize, for any e > 0, if Ro{e) > is such that |-Bro| = (1 — e) exp y — ^zii-j, then 

Mk ^ sRl 

k 

Since |Bro I = ^^ir^^o we get our conclusion with C{n, k,£)=£ (^h^) " • □ 

Remark A. 1. In dimension n = 3 with k = 2 and e = 1/2, one sees that there is some 
constant C > such that 

J^J-{v)\v\^dv ^ \ '''exp(-4H(/)/3) (A.4) 

for any nonnegative distribution function f(v) > with /j^s f{v)dv = 1. 
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